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Abstract 

A nonpolynomial one-dimensional quantum potential representing an oscillator, 
that can be considered as placed in the middle between the harmonic oscillator and 
the isotonic oscillator (harmonic oscillator with a centripetal barrier), is studied. 
First the general case, that depends of a parameter a, is considered and then a 
particular case is studied with great detail. It is proven that it is Schrodinger 
solvable and then the wave functions and the energies of the bound states 
are explicitly obtained. Finally it is proven that the solutions determine a family of 
orthogonal polynomials Vn{x) related with the Hermite polynomials and such that: 
(i) Every Vn is a linear combination of three Hermite polynomials, and (ii) They 
are orthogonal with respect to a new measure obtained by modifying the classic 
Hermite measure. 
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1 Introduction 



It is well known that, in quantum mechanics, the family of Schrodinger solvable potentials is very 
restricted and also that exact solvability is a very delicate property (see Ref. [T] for a review and 
Refs. p]-[llj for some recent papers on this matter). In fact, in most of cases, the addition of 
a small perturbation to a quantum solvable system breaks this property and leads to potentials 
that must be analyzed by the use of perturbative methods, variational formalisms or numerical 
techniques. The most interesting and best known system inside this small family is the harmonic 
oscillator whose energy spectrum consists of an infinite set of equidistant energy levels. Many 
other oscillators, as for example harmonic oscillators perturbed by a term containing a fourth 
or a sixth power in the coordinate, have been extensively studied but making use of the above 
mentioned techniques. Nevertheless, it is known the existence of another oscillator, the so-called 
isotonic oscillator, that is exactly solvable and is endowed with many interesting properties. 

The main aim of this paper is to present an study of a new nonpolynomial one-dimensional 
quantum potential representing an oscillator, that can be considered as placed in the middle 
between the harmonic oscillator and the isotonic oscillator. We will prove that is exactly solvable 
and that the energy spectrum and the wave functions have properties closely related with those 
characterising the harmonic oscillator. 

In more detail, the plan of the article is as follows: In Sec. II we recall the main characteristics 
of the isotonic oscillator In Sec. Ill we present a new potential depending of a parameter a; first 
we study the general case and then we solve the particular case = 1/2 obtaining the wave 
functions and energy spectrum. Sec. IV is devoted to analyze a family of orthogonal polynomials 
and to study the relation with the Hermite polynomials. Finally, in Sec. V we make some final 
comments. 



representing an harmonic oscillator with a centripetal barrier, is known as the isotonic oscillator 
\yi\ I13j . It is important because is endowed with properties closely related with those of the 
harmonic oscillator. At the classical level it leads to periodic solutions with the same period 
(isochronous potential [13]) and at the quantum level it is Schrodinger solvable, the Hamilto- 
nian is factorizable |15) and the energy spectrum is equidistant. Moreover, it is related with 
supersymmetric quantum mechanics |16j . The two-dimensional version is superintegrable and 
corresponds to the so-called Smorodinsky-Winternitz system and the centripetal term relates it 
with the Calogero-Moser system. In addition to all these theoretical properties, this particular 
nonlinear oscillator is important in quantum optics and in the theory of coherent states [171 118j . 



2 The Isotonic oscillator 



The following potential 




5>0 
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The classical Lagrange equation, that is given by 

x + w^x — ^ = 0, 

is a particular case of the Pinney-Ermakov equation [19]. It can be exactly solved and the 
solution is given by 
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X = — ^ y' {'^'^^'^ - g) sm^iyJt + <p) + g . 

showing explicitly the above mentioned periodicity. At the quantum level it is convenient to 
write g in the way g = m{m + 1) (the constant g should be greater than —1/4 and m can be 
any real number but we will always take it as non-negative) so that the Schrodinger equation 
takes the form 



T 

dx'^ 



2 2 , 
LO X + 



m{m + 1) 



x^ 



^ + 2E'^ = 0. 



where we assume h = 1 for easy of notation. The simplest solution, representing the ground 
state, is 

^'o = A^ox^+™ exp(-^wx2), EQ = {^ + m)uj. 
and the all the other wave functions and eigenenergies are given by 

^'2n = iV2nX^+™P2n(x) exp{- ]^UJ X^) , ^Sn = ^0 + 2n(J, 

where P2n{x) is a polynomial of order 2n with only even powers of x and A^2n denotes the 
normalization constant. The energy spectrum is equidistant since 

-E'2n+2 = E2n + 2 w , n = 0, 1, 2, . . . 

Nevertheless, the height IS.E of the energy steps is twice that of the simple harmonic oscillator 
Uq. In fact, it seems as if half of the levels (those with an odd number of nodes) have disappeared. 



3 A new solvable potential 



We now turn our attention to the study of the one-dimensional quantum system described by 
the following potential 



Uoa{x) = Uo{x) + Uaix) = {h 



u? x^ + 2ga 



x^ — 
'x^ + a^)^ 



5a > 0, 



where a is a positive real parameter. The reason for this particular algebraic expression is that 
the new additional term can be written as the sum of two centripetal barriers in the complex 
plane 



x2 - a2 1 1 

+ 



(x^ + a^)^ (x + ia)^ (x — ia)^ 
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so it is a rational potential with two imaginary poles symmetric with respect the origin. Actually, 
if ga remains constant, then when a goes to zero and to oo the potential Uoaix) converges to the 
isotonic and the harmonic oscillator, respectively. 

The first important property is that if the coefficient ga is not arbitrary but given by ga = 
2uja^{l + 2(^0^) then the Schrodinger equation, that takes the form 



■2 



^' T 



cj2x2 + 4a;a2(l + 2a;a2)-4 ^1^ + 2E^' = 0, (1) 



admits the following solution 

*o = p^f^ exp(-ia;x2), = \u: - {2u: af . (2) 

It represents the ground level (^'o has no nodes) and it clearly shows that when a ^ then 
the corresponding wave function and ground level energy of the linear oscillator is obtained. 
We also note that the above expression ga = 2a;a^(l + 2u;a^) can be alternatively written as 
9a = ruaima + 1) with nia = 2a; a^. 

In the following we will focus our study on the particular case = 1/2, nia = io (see Figure 
1) represented by the equation 



•2 



<^ T 



' 1 2x^-1 



(2x2 + 1)2 



^ + 2E^ = (3) 



so that ^0 and Eq become 



$o(x) = ^° ^ exp(- -x2) ^ Eo = --, (4) 
^ ^ 1 + 2x2 2 ^ 2 ' ^ ^ 

where we have assumed uj = 1 for easy of notation. 

As a method for obtaining all the other eigenstates we assume that the functions ^(x) can be 
factorized in the form 

^'(x) = F{x) ^'o(x) 

and, as the lowest energy is Eq = —3/2, it seems appropriate to write the general energy E as 
follows 

3 

E = -- + e. 

Then the Schrodinger equation ^ leads to the following equation for the function F{x) 

(1 + 2x2) ^// _2x{5 + 2x2) p' ^2e{l + 2x2) jp = q (5) 

Since the origin x = is an ordinary point we expect an analytic solution F with a power 
series expansion convergent in the interval {—R, R) with the radius of convergence R given by 
R = l/V2 



F{x) = ^ PriX" = Po+ PlX + P2X'^ + P3X^ + 



n=0 
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that when replaced in ([5]) leads to 



2(P2 + epo) + [6p3 + 2(e - 5)pi] x + ^ 



m=0 



m + 4)(m + 3)pm+4 

= . (6) 



+2[(m + 2)(m-4) +e]pm+2 + 4(e-m)p^ 

Therefore, we first obtain the following two relations for the first coefficients po,P2 and pi,P3 

P2 + epo = , 3p3 + (e - 5) pi = , 
and then the general recursion relation 

{m + 4:){m + 3)pm+4 + '2[{m + 2){m-i) + e]pm+2+ 4(e-m)p„ = 0, m = 0,l,2,... . 



Hence the recurrence relation involve three different terms {pm+A depends of both Pm+2 and Pr, 
and is constrained by the two first relations that are not included in the general rule. In any 
case even coefficients are related among themselves and the same is true for odd coefficients. 
The general solution will be obtained by fixing the values of po and pi, i.e. the value F{0) and 
F'{0). The solution determined by F{0) = 1, F'{0) = is even (only contains even powers) 
while the one determined by F{0) = 0, F'{0) = 1 only contains odd powers of x. Actually the 
expressions of the first coefficients, in terms of po and pi, are given by 

P2 = -epo = 2^(e _ 50)(e - 6)(e - 4)epo 
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P4 = jjie - W)epo pw = -^(e - 82)(e - 8)(e - 6)(e - 4)epo 
P6 = -i^(e-26)(e-4)epo 



and 



P3 = -jj,{e - 5)pi p9 = |l(e - 65)(e - 7)(e - 5)(e - 3)pi 

P5 = f,{e - 17)(e - 3)pi = -^(e - 101)(e - 9)(e - 7)(e - 5)(e - 3)pi 

P7 = -frCe - 37)(e - 5)(e - 3)pi 



In the particular case in which e is an even integer number e = 2k, with = 3, 4, 5, . . ., all the 
coefficients P2(k+r) vanish and the series reduces to a polynomial of degree k in powers of x"^, 
P2k{x)- Similarly, when e is an odd integer number e = 2A; + 1, with k = 2,3,4,..., all the 
coefficients P2{fc+r)+i with r > 1 vanish, and the series reduces to a polynomial P2k+i with only 
odd powers of x. The first polynomial solutions are: 

„2 A^A ,83 ^5 — X ^X 



=1-6X2-4X^ + |X'- _ 4^5, 8^7 

_ 1 _ e^2 _ 8 4 I 32 6 _ 16 8 /^r - X ^X ^X +21^ 
-Ts — -L OX gX -h 5 X ^gX _ „ 4 3 8 ^5 , 16 ^7 16 „9 

PlO = 1 - 10x2 + fx6 - |?X« + -^xlO -^-3^-5^+15^-135^ 



10 = 1 - lux- -I- -g-x- - 2iX- -I- 

Two important properties are: First, there are no polynomial solutions of degree k = 1 and 
k = 2 and, second, all of them have two complex conjugate roots, so that P3 has an unique real 
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zero, P4 has two real zeros and, in the general case, P2fc and P2k+i with k > 1 have 2k — 2 and 
2k — 1 real zeros respectively. 

Another remarkable property is that when the polynomials Pn are expressed as linear combi- 
nation of Hermite polynomials then only a finite number of terms are different from zero. We 
have obtained the following relations for the first cases 



P4 = -4Ho - 4.H2 - P3 =Hi + ^Hs 

Pq = 3H2 + H4 + P5 = —Hi — — 5)-f^5 

Ps = -IH4 - ^Hq - P7 = + + gigi^T 

Pio = j^He + ^Hs + ^Hio P9 = -^Hr, - ^^Hj - 

These particular relations clearly suggest that each polynomial P2k can be written as a linear 
combination of H2k and the two previous even Hermite polynomials, H2k-2 and i?2fc-4; similarly 
the odd polynomial P2k+i appears as a linear combination of only H2k+i, -f^2fc-i and i?2fc-3- 



4 Wave functions and orthogonality relations 

The differential equation dS]) is not in selfadjoint form because, if we denote by ao, ai and 02 the 
three coefficients of the equation, we have 

ao = 1 + ai = -2x (5 + 2x'^) , / ai . 

ax 

However, use can be made of an integrating factor fi{x) such that 

d 



so that is given by 



— [fi{x)ao{x)] = fi{x)ai{x) . 
ax 



oo (1 + 2 x^j'^ 



in such a way that ([5]) becomes 



dx 



p{x[ 



dF^ 

dx 



+ 2er{x)F = 0, 



where the two functions p = p{x) and r = r{x) are given by 



p(x) 



(1 + 2x2)2 ' 



r{x) 



(1 + 2x2)2 



This equation, with the appropriate conditions for the behaviour of the solutions at the end 
points, constitutes a Sturm-Liouville (S-L) problem, defined in the real line M. According to 
this, the eigenfunctions of the S-L problem are orthogonal with respect to the weight function 
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r = e ^/(l + 2x^)^, and, in particular, the polynomial solutions Pm, m 
differential equation ([5]), satisfy the orthogonality conditions 



0,3,4,..., of the 



(1 + 2x2)2 



= , m ^ n . 



Let us now return to the Eq. ([5]) and suppose for F the following factorization 

F' = {l + 2x^)G. 

Then we arrive after some calculus (we omit the details) to the following equation for the function 
G 

G" - 2x G' + 2(e - 3) G = , 
that means that the derivative -P^(x) of the polynomial P„(x) must satisfy 

= (1 + 2x2) /7„_3 , e = n, n = 3,4,5,... 

(up to a multiplicative constant). At this point we recall the following two properties of the 
Hermite polynomials 



ii) 



2xH„ 



Hm+i + 2mH, 



m—l 



{a) H'^ = 2mHm-i 
Then making use of (i) we arrive to 



P^, = ^ [H,,-i + 4(n - 2) /f„_3 + 4(n - 3)(n - 4) H^^^] , 



and making use of (ii) and then integrating we obtain 

1 



An 



Hn + (4n) Hn-2 + (4n)(n - 3) Hn-i 



It seems convenient to multiply P„ by 4n and introduce the new family of polynomials Vn defined 
in the form 

Vn = Hn + AnHn^2 + ^n{n-?,)Hn-i, n = 3,4,5,... (7) 
so that the coefficient of Hn (the dominant term in the expression of Vn) reduces to unity. 



Proposition 1 The following equality holds 



2nPne-^ 

(1 + 2x2)2 



d_ 

dx 



Hn-3 ^-x2 



1 + 2x2 



, n = 3, 4, 5, . . . 
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Proof: This statement is proven just by making the calculus. 
Now we can write 



J —I 



(1 + 2x2)2 

and integrating by parts wc arrive to 



1 d 

(16nm) I — 

^ ^ '2mdx 



m—3 



1 + 2x2 



Pn{x) dx 



I 

J —(. 



(1 + 2x 



2^2 



dx 



6n 



m—3 — 

e 



1 + 2x2 



J — oo ^ I ^■^ 

/oo 
-oo 



PLix) dx 



3 6 ^ dx 



= <5^„(8n)[2--=^(n-3)!v^] 



(8) 



So, we relate the problem of normalization with the standard problem of Hermite polynomials. 
Hence the orthogonality conditions for the family Vn{x) read: 

/ Vm{x) Vnix) r(x) dx = d^n / . ^0 U dx = kn {2^ n\ V^) , 

J-oo J -oo + ^X ) 

where the proportionality constant kn is given by 

1 

kn 



(n- l)(n-2) ■ 

Sec Figures 2a and 2b for the plot of some of these polynomials. 
If we define the P-Hermite functions by 



'Pmjx) (1/2)^2 

(1 + 2x2) 



, m = 0, 3, 4, . . . 



then the above property admits the following alternative form: 

"^rnix) ^^nix) dx = , 171 ^ n . 



In summary, the eigenfunctions corresponding to the lowest energy levels are: 

Voix) 



*o(x) 


= No 


*3(a;) 




■^4{x) 


= N4 


^5{X) 





(1 + 2x2) 
Vsix) 

(1 + 2x2) 
Vi{x) 

(1 + 2x2) 

n(^) 

(1 + 2x2) 



e-(V2)-% £;o = -3/2 

£;3 = 3/2 = -3/2 + 3 
g-(i/2)x2^ £;4 = 5/2 = -3/2 + 4 
g-(i/2)x2^ £;5 = 7/2 = -3/2 + 5 
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where the normahzation constant is 



Nk 



{k-l){k 



1/2 



The first three wave functions, ^o{x), ^^(x), ^4(x), together with the corresponding wave 
functions of the harmonic oscillator ^o{x), ^i{x), ^2ix), are plotted in Figures 3, 4 and 5. 

The energy £"0 of the ground state ^o{x) has been singled out of all the other values and 
moved into the smaller value Eq = — 3/2. The rest of the energy spectrum consists, as in the 
pure harmonic case, of an infinite set of equidistant energy levels 



n+l 



En + l 



n 



3,4,5,... 



Let us close this section with two comments on the new family of polynomials we have obtained. 
First the definition ([7]) of Vn as a linear combination of only three Hermite polynomials can be 
considered as a particular case of a situation known as a special linear combinations of orthogonal 
polynomials (see [201121) and references therein). Finally, let us mention that taking into account 
the "Rodrigues formula" for the Hermite polynomials 



-1 



dx'^' 



we obtain 



rn{x) = i-iy 



+ 4n , ^ ^ + 4 n (n - 3) 



d 



'ra— 4 



dx^ dx^~'^ dx"-~^ 
that must be considered as the "Rodrigues formula" for this new family of orthogonal polyno- 
mials. 



5 Final comments and outlook 

We have proved that the potential U^aix) can be exactly solved in the particular case a? = 1/2 
and also that it possesses two very remarkable properties. First, the fundamental level ^'o has 
an energy Eq that is lower than in the pure harmonic case and, in a sense, is isolated of all 
the other values. Second, the rest of the energy spectrum is endowed with the equidistance 
property. Concerning the general case, with an arbitrary value for the parameter a, we have 
only obtained the expression for the fundamental level ^{a) , EQ[a)) . The resolution of the 
general case remains as an open question that deserves be studied. Finally let us also mention 
that the analysis of this potential using the supersymmetric quantum mechanics as an approach 
also seems an interesting matter to be studied. 
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Figure Captions 



• Figure 1. Plot of the potential Uoa{x) for cj = 1 as a function of x for = 1/2 
(continuous line) together with the plot of the harmonic oscillator (dash line). The main 
difference lies in the form of the minimum that is much deeper in the Uoaix) case than 
in the linear case. Nevertheless, for great values of |x| the two functions have rather the 
same form. 

• Figure 2. (2a) Polynomials V3 (dash line) and V4 (continuous line). P3 has a unique 

real zero at the origin and P4 has two real zeros (symmetric with respect the origin). (2b) 
Polynomials V5 (dash line) and Vq (continuous line). P5 has three real zeros (the origin 
and two other placed symmetric) and Pg has four zeros (two positive and two negative).. 

• Figure 3. Wave function ^'o (continuous line) and wave function $0 of the harmonic 

oscillator (dash line). 

• Figure 4. Wave function ^'3 (continuous line) and wave function $1 of the harmonic 

oscillator (dash line). 

• Figure 5. Wave function ^'4 (continuous line) and wave function $2 of the harmonic 
oscillator (dash line). 
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Figure 1: Plot of the potential Uoa{x) for a; = 1 as a function of x for = 1/2 (continuous 
line) together with the plot of the harmonic oscillator (dash line). The main difference 
lies in the form of the minimum that is much deeper in the Uoa{x) case than in the linear 
case. Nevertheless, for great values of |a;| the two functions have rather the same form. 




Figure 2: (2a) Polynomials V3 (dash hue) and V4 (continuous line). P3 has a unique real 
zero at the origin and P4 has two real zeros (symmetric with respect the origin). (2b) 
Polynomials (dash line) and Vq (continuous line). P5 has three real zeros (the origin 
and two other placed symmetric) and Pe has four zeros (two positive and two negative). 




Figure 3: Wave function ^0 (continuous line) and wave function $0 of the harmonic 
oscillator (dash hue). 
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Figure 4: Wave function ^3 (continuous line) and wave function $1 of the harmonic 
oscillator (dash line). 




Figure 5: Wave function ^4 (continuous line) and wave function $2 of the harmonic 
oscillator (dash line). 
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